Ab initio determination of the lifetime of the Q 2 Pz/2 state for 207 Pb + by relativistic 
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Relativistic coupled-cluster(RCC) theory has been employed to calculate the lifetime of the 6 2 P3/2 
state of sins le ionized lead( 2(17 P6) to an accuracy of 3% and compared with the corresponding 
value obtained using second order relativistic many-body perturbation theory (RMBPT). This is 
one of the very few applications of this theory to excited state properties of heavy atomic systems. 
Contributions from the different electron correlation effects are given explicitly. 
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I. INTRODUCTION 

Trapped and laser cooled ions are excellent candidates for 
high precision measurements PIQ- 207 Pb + is the heaviest 
atomic ion to be trapped and cooled to date0,0- The 
lifetime of the 6p3/2(6 2 P 3 / 2 ) state of this ion has been 
measured to an accuracy of 2%@. The transition from 
this excited state to the ground state 6pi/ 2 (6 2 Pi/ 2 ) is pre- 
dominantly magnetic dipole(Ml) in character, but there 
is a small electric quadrupole(E2) component as well. 
The Ml and E2 transition rates depend on the inverse 
cubic and quintic powers of the transition wavelength re- 
spectively. The lifetime of the 6 2 P3/ 2 for Pb + must there- 
fore be calculated by an accurate relativistic many-body 
method. Indeed a comparision of the measured and cal- 
culated values of the lifetime for this state would consti- 
tute a stringent test of the theoretical approach on which 
the calculation is based. In addition, a precise knowledge 
of this lifetime is useful in determiningthe abundance of 
Pb in the planetary nebula NGC 70270, Q - The work in 
this paper uses the non-linear RCC theory to calculate 
this lifetime. This theory has been successfully applied 
earlier to different atomic systems 0, 0, 0, [n|, but it 
has seldom been used to study excited state properties 
for heavy atoms or ions. This is the first accurate calcu- 
lation of the lifetime of the 6p 2 P 3 / 2 state of Pb + to our 
knowledge. We compare our RCC results for the excited 
energy and lifetime with those obtained from second or- 
der RMBPT(2). 
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(2.1) 



Uf\\Q^\\ji) = (jf\\C^\ji)^ 
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respectively, where ji and ft, are the total angular 
momentum and relativistic angular momcntum(K, = 
±(j'i + |)) quantum numbers respectively of the i'th 
electron orbital and k = Lua, where lu — e/ — q, is the 
photon energy of the transition, a is the fine-structure 
constant. We use atomic units (% — m e = |e| = 1) in 
this paper. The quantity Cm\f) is the Racah tensor 
and ji(kr) is the spherical Bessel function of order 
I. P Ki and Q Ki are the large and small radial compo- 
nents of the Dirac-Fock i th single particle wave functions. 

The angular factor is given by 



< jf\\C^\\ji >= (-1)^ +1/2 Vj/ + 1/2v^TT72 
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II. THEORY 

The one-electron reduced matrix elements of Ml and E2 
operators are given by 01 

OLK I 



The Ml and E2 transition probabilities Api(s~ 1 ) can be 
expressed in terms of the line strength Sfi which is the 
square of the sum of the single particle transition matrix 
elements given by eq(2.1) and eq(2.2) for the appropriate 
transitions and wavelength \(A°) as [l3l Il4| 
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and 



E2 _ 1.1199 x 10 18 QE2 
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where J/ is the degeneracy of the initial metastable state 
which is equal to 3/2 for the present calculation. 



III. METHOD OF CALCULATION : 
RELATIVISTIC COUPLED CLUSTER THEORY 

The relativistic Dirac- Coulomb atomic Hamiltonian is 
given by 

H = V ca.fi, + (/3- l)c 2 + V nuc ( rj ) + V — (3.1) 



j 



3<l 



where a and (3 are the usual Dirac matrices and V nuc (rj) 
is the potential at the site of the j th electron due to the 
atomic nucleus. The energy eigen values are scaled with 
respect to the rest mass energy of the electron. We first 
solve the relativistic Hartree-Fock(Dirac-Fock) equations 
to obtain the single particle orbitals and their energies. 

H D f = cd -PJ + s - l )° 2 + Vnuc(rj) + U DF (rj) 



The residual Coulomb interaction is given by 
V. 



3<l J 3 



The single particle orbitals are obtained by solving the 
following equation self-consistently 



where 



and 



{t j + U DF (r i ))\4> i )=e j \4> i ) 



tj = ca ■ pj + {(3- l)c 2 + V nuc (r 3 ) 



(3.3) 
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The single particle relativistic orbitals can be expressed 
as 

U.( r )\ -If P s< r )\x Kjmj ) \ 

where Pj(r) and Qj(r) are the radial part of the large 
and small components respectively and |x«<m<) and 
\X-K-m ) ar e their respective spin angular momentum 
components. Cj's are the single particle energies. 

We have employed the RCC to incorporate correlation ef- 
fects among electrons due to the residual Coulomb inter- 
action. In this approach the exact atomic wavefunction 
for the closed-shell system can be expressed as[l7| 



|*cc) = e T |*) 



(3.4) 



where T is the core electron excitation operator and |$) 
is the above closed-shell dctcrminantal state built out of 
the Dirac-Fock single particle orbitals. 

In the closed-shell coupled-cluster theory one starts with 
the equation 

He T \^) = Ee T \<S>) (3.5) 

The energy and amplitude determining equations are 

=E5 Kfi (3.6) 

where H = e~ T He T , \& K ) is a determinantal state 
with K = 0,1,2.... representing the reference state and 
excited determinantal states. We have considered all 
possible non-linear terms in T- operator for its amplitude 
determining equations. 



Goldstone[l7l [ l^ . Il9l | and angular momentum 
diagrammatic [I3 |2(| techniques are used for evalu- 
ating different radial integrals and angular factors. The 
normal ordered Hamiltonian is defined as 



H N = H - (*\H\$) = H -E DF , (3.7) 
where E DF = 

We have truncated our wavefunction expansion at the 
level of singles and doubles(CCSD) and all possible non- 
linear terms have been included in the above equation. 

The ground state of Pb + contains only one valence 
electron; namely the 6pi/2 orbital. One way to evaluate 
the ground state energy of Pb + is to first compute the 
wavefunctions for the closed shell system Pb ++ using 
the above closed shell CC approach and then append a 
valence electron (6pi/ 2 orbital) using the open shell CC 
(OSCC) method as follows. The same procedure has 
been followed to obtain the excited 6p 3 / 2 state. 

The new reference state of the open-shell system with 
one valence electron v can be expressed as|21| 

\$ v ) = 41$) (3.8) 

where al is the particle creation operator. The ex- 
act atomic states are defined now, using the Fock-space 
OSCC method, as [Hip 



|^)=e T {e s "}|$„) 



(3.9) 



where S v is the valence excitation operator. Since the 
system under consideration has only one valence electron, 
the S- operator exponetial series naturaly truncates at 
linear term, i.e. the open-shell wavefunction has the form 



|tf„) =e T {l + S v }\f> v ) 



(3.10) 



where 



S v = S vl + S v2 = y^^ Qp ays? + 77 ^ at a t a b a v s v. 
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and 



S v i — S a^a v s^ 



Sv2 = \Yl a : 



p a q CLbCL v S vb 



(3.11) 



bpq 



with s"P and are the cluster amplitudes corresponding 
to single and double excitations involving the valence 
electron. 

In the next step, we include approximate triple exci- 
tations by contracting the two-body operator's) and 
the double excitation operators (T 2 , 5^2 ) in the following 
wayHIH 



qpqr _ 
°vbc — 



V es T2 + V es S v 2 



(3.12) 



where e$ is the orbital energy of the i'th orbital. Note 
that we use notations a,b,c..., p,q,... and i,j,... for the 
core, particle and generic orbitals respectively. 

The wavefunction in the framework of many-body per- 
turbation theory(MBPT) can be written as 



I* 



AIBPT 



} = + + |$ (2) ) + |$ (3) ) + ... (3.13) 



where | is the n th order correction to the wavefunc- 
tion |$). Each order corrected wavefunction is a linear 
combination of all excited determinantal states with re- 
spect to |$). Gathering excitations of the same order 
together from each of the corrected wavefunction we can 
rewrite the wavefunction of the many-body system as 

\*mbpt) = + X;(C*(1, 1) + C*(2, 1) + ..)|*S) 

ap 

= |$„> + X;(CJ(1 ) 1) + CJ(2,1) + ..)I*5> 



p^u 



+^£( C £( 1 ' 2 ) + C S6( 2 ' 2 H 

abpq 
bpq 

\*v)+ £(CS(1,1) +CS(2,1)- 
+ £(C*(1,1) + C?(2,1) 



■)I*S> 



l*S> 



.)<4°a|*> 

.)ala„|$) 



p^u 



i £)(C™(1, 2) + C™(2, 2) + ..)at a t flbaQ |$) 



abpq 



\ ^(^(1,2) + Cft(2, 2) + ..) a t a t abat ,|^) 

6pg 



bpq 

+ \ £ tlt\a\a\a h a a $) + 



abpq 



1 



(1 + Ti + T 2 + ~I?...)|*„) + Si„|*„) + S 2 „|$„) 



= e T {l + ^}|<i> !J } 
s |*cc) 



(3.14) 



where 



£ = C*(1,1) + C*(2,1) + .... 
< = C£(1,1) + C£(2,1) + .... 

C = ^(l,2) + C^(2,2) + .... 

Tl = *a a p a « 
ap 

Si„ = ^ s£a p a a 
p^d 

&. = *£« 



H lapa* q a b a v 



bpq 



and 



T = Ti + T 2 

ap abpq 

S V = Sly + 02 V ' 

^ s^apa^ + - Y s p v q b ala\a b a v 



(3.15) 



p^u 



bpq 



and C(l, 1), C(l, 2), ...etc. are the perturbation co- 
efficients for each order of the corrected wavefunction. 
The first index of the superscript represents the order of 
the pertubation and the second represents the excitation 
level. Therefore, each T- operator accounts for the 
correlation effects from all orders of the perturbed 
wavefuction. The above relation shows that unlike in the 
MBPT, where the correction to correlation is computed 
order by order, in coupled-cluster theory once the 
leading correlation corrections are identified, a subset 
of terms to all orders which improve these corrections 
are also included. Diagrammatic representation for this 
theory has been shown in the Fig. 1. 

The equations for the open-shell cluster amplitudes are 
determined from |l7ll21| 



($ V \H N {1 + S v }\$ v ) = AE{v) 



(3.16) 



and 



abpq 



(K\H~n{1 + S v }\$ v ) = —AE(v)($*\{S v }\$ v ) (3.17) 



4 



A P 



P^ v 

lv 





A 





2v 

FIG. 1: Factorisation of all order CC amplitude diagrams for 5i„ and Sbu operators in terms of counterpart MBPT diagrams 



where, AE(v) is the electon attachment energy which is 
equal to the negative of the ionisation potential for the 
valence electron, v. 



If 



and 



|*,) =e T {l + S v }\$ v ) 



|1V) = e T {l + 5„/}|$„») 



represent the ground and excited states with energies 
E(v) and E{v') respectively then the excitation energy 
is given by 



E(v') - E(v) = AE{v') - AE(v) 



(3.18) 



The transition matrix element for a general one particle 
operator can be expressed in coupled-cluster theory as 



Ofi = 



<*/|0|*i> 



(*f\0\%) 



y/(l + N f )(l+Ni) 

+ S}}e Tt Oe T {l + S z }\$ l ) 

y/{l+N f )(l + Ni) 



(3.19) 



+ 5}}Q{1 + S i }|$ i ) 
V /(1 + A r /)(1 + ^) 



(3.20) 



where the normalisation terms for the V th orbital is ob- 
tained from 

N v = ($ v \Si[e Ti e T ] + Sl[e T 'e T ]S v + [e Tt e T ]^|$„) 
= ($ v \Sifi v + Sln v S v + n v Si\$ v ) (3.21) 



with 



O = (e 1 Oe') f . c . + {(■/ ()c/) . b . + (e T+ Oe r ) t . b . + 



(3.22) 



and 



n v = {e T e T ) f . c . + (e T e T ) . b . + (e T ' e T ) t . b . + ... (3.23) 

The f.c, o.b., t.b.,..etc abbreviations are used for the 
fully contracted, effective one-body, effective two-body 
...etc terms respectively [24|. Terms containing only upto 
effective three-body diagrams will contribute to both the 
numerator and the denominator. The fully contracted 
terms are excluded on the basis of the linked-Diagram 
theorem^?] in the evaluation of the O and N. All the 
one-body terms have been taken into account as their 
contribution to the correlation effects is the largest. The 
dominat par ts of the two-body terms have also been 
computed|lOl |24| . Finally, these terms are contracted 
with S f and Si operators. 
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Contributions from the normalisation factor have been 
determined in the following way 



Norm = (*/|0|*i){ 



1 



y/il+Nf^l+Ni) 



- 1} (3.24) 



IV. RESULTS AND DISCISSIONS: 

We have used Gaussian type orbitals(GTO) for the con- 
struction of single paricle orbitals of the Dirac-Fock 
wavefunction(|$}), whose expression is given bv|l5| 



matrix elements, the total contribution of OSu and its 
conjugate is larger than OS21 and its conjugate term. In 
its lowest order OSu corresponds to the Brueckner pair 
correlation and OS^ to core polarisation. The largest 
contribution to electron correlation comes from the 
pair correlation effects. The reduced transition matrix 
elements for Ml and E2 operators are given in table II. 
The important contributions from the one-body terms 
of these two quantities are given in table III. 

The net probability for a given transition which allows 
two different channels is given by 01 



F; 



(L/S) 



(r) 



E(L/S) k - 
c\ >r e 



(4.1) 



A = A l 
This can be expressed as 



with k=0,l,2,... for s, p, d,- ■ ■ respectively. The function 
F^ S \r) stands for the large(L) and small(S) compo- 



1 

n 



A 2 
1 

T 2 



(4.3) 



(4.4) 



nents of the dirac wavefunction. c\ L ^ S ^ is the expansion 
coefficient of the corresponding large and small compo- 
nents respectively. The kinetic balance condition |25j has 
been imposed between the large and small components 
of the GTOs. For the exponents, the even tempering 
condition 



i/3, 



1. 



■,N (4.2) 



has been applied. Here, N stands for the total number 
of basis functions for a specific symmetry. In the present 
calculation, we have taken a = 0.00825 and (3 — 2.73 
for all symmetries. We have considered 13si/ 2 , 13pi/ 2 , 



where T\ and r 2 are the lifetimes through different 
branches and r is the total lifetime. The results obtained 
for Vj which represents the Ml lifetime and 't' 2 the E2 
lifetime are given in table IV. The lifetime for the 6 2 Pj/ 2 
state obtained using RCC is 0.0425 second which has an 
accuracy about 3%. The lifetime result of our second or- 
der RMBPT calculation is 0.0418 second. This suggests 
that there would be a strong cancellation of the correla- 
tion effects from the higher order RCC contributions. 



13P; 



'3/2, 



lid 



3/2, 



lid 



5/2, 



15/2: 



8/7/2, 7 97/2 and 7g g / 2 
active orbitals. All core electrons have been excited in 
the present calculation. 



V. CONCLUSION 

Relativistic coupled-cluster theory has been applied to 
calculate ionisation potential and transition matrix ele- 
ments for Ml and E2 operators to determine the lifetime 
We have obtained an accuracy of better than one percent for the fi2p 3/ 2 state of the singly ionised lead. The accu- 
for the ionisation potentials of both 6 2 P 1/2 and 6 2 P 3/2 rac y of the calculation is about 3%. 
states. The ionisation potential at the Dirac-Fock level 
is calculated using Koopman's theorem. As can be seen 
from table I the correlation contributions are about 5% 
and 6% for the 6 2 Pi/ 2 and the 6 2 P3/ 2 states respectively. 
However, the CC excitation energy improves only by 
about one and half percent over the DF value and 
deviates by 2.7% from the experimental value. The 
excitation energy between these two states is calculated 
using eq. (3.18) and it's value is given in table IV. From 
table II, it is clear that for both the Ml and E2 transition 
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TABLE I: Ionization potential energies of different states of Pb + 




States 


Koopman(cm- 1 ) MBPT(2)(cm- i ) CCSD(T)(cm- i ) Expt^cm" 1 ) 


% of accuracy 


6f>l/2 

6P3/2 


114015 121898 120126 121208 
100402 108041 106416 107123 


0.8 
0.6 






TABLE II: Contributions from the important terms to the Ml and E2 transition matrix elements in a.u.. 




(i) Contribution from one-body terms: 






Terms 


Contributions for Ml 


Contributions for E2 


6 

OSu 

s\ f o 
0s 2i 
stp 

S 1 ^OSu 

SlfOS M 
Norm. 


-1.12021925042 
-0.00517412601 
0.00553162299 
0.00017425249 
-0.00031546156 
-0.00112434477 
-0.02963327056 
0.00989789079 


9.05230735868 
-0.18227019838 
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-0.09056241322 
0.08833942426 
0.01932259992 
0.19544605525 
-.07635041195 



(ii) Contribution from two-body terms: 



Terms Contributions for Ml Contributions for E2 
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Wavelength 
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E2- Channel 


Total 




(cm" 1 ) 


(in A) 


(n) 


(r 2 ) 


(t) 


Dirac-Fock 


13613 


7346 


0.0455 


0.909 


0.0433 


MBPT(2) 
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7216 


0.0436 


1.03 


0.0418 


CCSD(T) 


13710 


7294 


0.0445 


0.958 


0.0425 


Exptf5l 


14085 


7100 






0.0412(7) 
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